Interdisciplinary Group

• Computer Science
• Applied Statistics and Scientific Computation
• Computational Biology

Central theme

• Computational and statistical method and tool development
• Study the molecular basis of variation in development and disease
• Using high-throughput experimental methods

Metagenomics

• DNA sequencing of bacterial communities
• Method and tool development: metagenomeSeq Bioconductor package
• Large epidemiological study of childhood diahrrea in developing countries

Host-pathogen interaction

• Joint RNA-sequencing of parasite and host cell during course of infection
• Novel methodology for RNA-seq normalization and differential expression analysis

Computational Epigenomics

NHGRI strategic plan

NHGRI strategic plan

"The major bottleneck in genome sequencing is no longer data generation—the computational challenges around data analysis, display and integration are now rate limiting. New approaches and methods are required to meet these challenges."

• Data analysis
• Visualization
• Data integration
• Computational tools and infrastructure

Genes are expressed differently during different stages and in different tissues.

DNA is packed, making certain parts inaccessible, and this packing is dynamic.

DNA methylation is a chemical modification of DNA, involved in gene expression regulation.

Large blocks of hypo-methylation in colon cancer

• overlaps with other important genomic domains
• tissue-specific genes are over-represented within blocks

Hypo-methylation blocks observed across five solid tumor types.

Methylation pattern reconstruction problem

methylFlow: http://github.com/hcorrada/methylFlow

Methylation pattern reconstruction problem

• Given a set of mapped reads

Methylation pattern reconstruction problem

• Given a set of mapped reads
• Determine composition of cell-specific methylation patterns

Methylation pattern reconstruction problem

The statistic: number of reads in genomic region

• regions (vertices) defined by overlapping reads with consistent methylation patterns
• overlaps (edges) defined by overlapping reads with inconsistent methylation patterns
• region coverage: total number of reads originating in region

The model: expected number of reads in genomic region

\[ \mathbb{E} y_v = \sum_{u:(v,u) |in E} \ell_{vu} \sum_{p:(v,u)\in p} \theta_p \]

The estimator

• Penalized method of moments:
• number of parameters = number of paths through graph
• sparsity inducing penalty to obtain solution with small number of patterns

\[ \min_{\theta_p} \sum_v \lvert y_v - \sum_{u:(v,u)\in E} \ell_{vu} \sum_{p:(v,u)\in p} \theta_p \rvert + \lambda \sum_p \lvert \theta_p \rvert \]

How to solve efficiently

\[ \min_{\theta_p} \sum_v \lvert y_v - \sum_{u:(v,u)\in E} \ell_{vu} \sum_{p:(v,u)\in p} \theta_p \rvert + \lambda \sum_p \lvert \theta_p \rvert \]

If we interpret abundance as path flow, then we can rewrite in terms of edge flows

\[ f_{vu} = \sum_{p:(v,u) \in p} \theta_p \]

How to solve efficiently

\[ \min_{\theta_p} \sum_v \lvert y_v - \sum_{u:(v,u)\in E} \ell_{vu} \sum_{p:(v,u)\in p} \theta_p \rvert + \lambda \sum_p \lvert \theta_p \rvert \]

If we interpret abundance as path flow, then we can rewrite in terms of edge flows

\[ f_{vu} = \sum_{p:(v,u) \in p} \theta_p \]

\[ \min_{f \geq 0} \sum_v \lvert y_v - \sum_{u:(v,u)\in E} \ell_{vu} f_{vu} \rvert + \lambda f_{vt} \\ \textrm{s.t} \sum_{u:(v,u) \in E} f_{vu} = \sum_{w:(w,v) \in E} f_{wv} \]

How to solve efficiently

\[ \min_{f \geq 0} \sum_v \lvert y_v - \sum_{u:(v,u)\in E} \ell_{vu} f_{vu} \rvert + \lambda f_{vt} \\ \textrm{s.t} \sum_{u:(v,u) \in E} f_{vu} = \sum_{w:(w,v) \in E} f_{wv} \]

• This is a linear optimization problem (LP)
• We can solve efficiently for very large problems
• Final solution obtained by path flow decomposition

Pattern reconstruction from whole genome bisulfite sequencing

Dataset of 50bp reads from mouse wild-type activated B cells, two types of progenitor cells (CLP and KSL).

Reconstruct patterns 4-100x basepair length

Pattern reconstruction from whole genome bisulfite sequencing

Reconstruct patterns with accurate marginal estimates

Pattern reconstruction from targeted bisulfite sequencing

Compare patterns across samples and populations

• We are now able to estimate methylation pattern composition for a single sample (e.g., normal or tumor)
• How to detect differences between cell populations:
  • For paired normal-tumor data, find genomic regions where methylation pattern composition changes (significantly)?
  • For samples across developmental course, how does composition change as differentiation occurs?
• How should we think about population-level inferences from measurements that are themselves inferred?

• New view into molecular profiling
  • Complex relationship between cell-to-cell differences and average methylation differences
  • Complementary to single-cell bisulfite sequencing
• Efficient formulation as network flow problem
  • Exploring relationship to stochastic optimization
• Code: https://github.com/hcorrada/methylFlow/
  • Includes incipient R package to parse and analyze resulting data
  • Takes SAM/BAM file input from Bismark and BSMap
  • Happy to collaborate on new analysis if you want to try it out
• Why do we care about intra-tumor heterogeneity?

Genes with hyper-variable expression in colon cancer are enriched within these blocks.

Gene expression hyper-variability enriched in hypo-methylation blocks in other cancer types.

Genes with consistent hyper-variable expression across tumors are tissue-specific.

• molecular methods for cancer detection, prognosis and treatment matching will be the basis of individualized medicine
• gene expression profile methods have been subject of study for decades
• very few proposed predictors are translated to the clinic
• one of the biggest culprits is lack of replicability of results in preliminary studies

anti-profile score: measures sample-specific deviation from normal expression in consistently hyper-variable genes

Anti-profile score distinguishes between stages in tumor progression

DNA methylation anti-profiles score distinguishes between stages in tumor progression

Stratification based on anti-profile score

Stratification of breast samples based on anti-profile score

• Profiles learned based on hyper-variability show consistent behavior across tissues and across experiments in tumor prognosis and progression
• Distinguish observations from two anomalous groups (e.g., adenoma vs. tumor)
• How can we incorporate the fact that we are classifying anomalies?
• Incorporating normal samples when building anomaly predictors improves stability and prediction performance
• Why (and when) is it worth doing that?
• Using function approximation methods to study predictor stability

• better understand connection between intra-tumor heterogeneity and consistent hyper-variability in cancer
• how to understand population-level inferences from features (measurements) that are themselves inferred
• move anti-profiles closer to the clinic
• explore anomaly classification as a general learning setting
• methods to understand hierarchical organization of epigenomic domains

• Discovery: consistent hypo-methylation and hyper-variability associated with progression and prognosis
• Methods: anomaly classification as a setting to understand predictor stability, methylation pattern reconstruction
• Tools

Collaborative and exploratory analysis

• Data transformation and modeling: data smoothing, region finding (R/Bioconductor: Bsmooth, minfi)
• Genome browsing: search by gene, search by overlap
• Region analysis: overlap with other data (our own, other labs, UCSC, ensembl)
• Regulation: expression data (Gene Expression Barcode)

• We have unprecedented ability to measure
• and lots of publicly available data to contextualize it

Integrative, visual and computational exploratory analysis of genomic data

• Browser-based
• Interactive
• Integration of data
• Reproducible dissemination
• Communication with R/Bioconductor: epivizr package

Creativity in exploration

We are building software systems to support creative exploratory analysis of large genome-wide datasets…

Computed Measurements: create new measurements from integrated measurements and visualize

Summarization: summarize integrated measurements (computed on data subsets)

Statistically-guided exploration: Calculate a statistic of interest

# Get tumor methylation base-pair data
m <- assay(se)[,"tumor"]

# Compute regions with highest variability across cpgs
region_stat <- calcWindowStat(m, step=25, window=80, stat=rowSds)
s <- region_stat[,"stat"]

Explore data based on statistic

What's around the regions with highest across CpG variability?

# get locations in decreasing order
o <- order(s, decreasing=TRUE)
indices <- region_stat[o, "indices"]
slideShowRegions <- rowRanges(se)[indices] + 1250000L
mgr$slideshow(slideShowRegions)

dynamically extensible: Easily integrate new data types and add new visualizations.

• Based on classic "three-table" design in genomic data analysis
• Data providers define coordinate space

One interpretation of Big Data is Many relevant sources of contextual data

• Easily access/integrate contextual data
• Driven by exploratory analysis of immediate data
• Iterative process
• Visual and computational exploration go hand in hand

Visualization goals

• Context
  • Integrate and align multiple data sources; navigate; search
  • Connect: brushing
  • Encode: map visualization properties to data on the fly
  • Reconfigure: multiple views of the same data

Visualization goals

• Data
  • Select and filter: tight-knit integration with R/Bioconductor;
  • (current work) filters on visualization propagate to data environment
• Model
  • New 'measurements' the result of modeling; perhaps suggested by data context

• collaborative computational and visual analysis
• effective visual methods to explore hierarchical organization of (epi)-genome
• deeper integration of statistically-informed visualization
• visualization-informed statistical analysis

• Discoveries: consistent hypo-methylation, hyper-variability
• Methods: anomaly classification as a setting to understand predictor stability
• Tools: computational and visual exploratory genomic data analysis

NHGRI strategic plan

"Meeting the computational challenges for genomics requires scientists with expertise in biology as well as in informatics, computer science, mathematics, statistics and/or engineering."

A new generation of investigators who are proficient in two or more of these fields must be trained and supported.

Courses taught and/or created:

• Undergraduate Computational Biology: sequence analysis and beyond
• Graduate Functional Genomics
• Undergraduate Data Science

Acknowledgements

Past members of HCBravo group
now at Harvard, U. Chicago, Johns Hopkins, Genentech, Dow Jones Data Science

Colleagues at CBCB
Current members of HCBravo group
Collaborators at JHU/Harvard

Funding: NIH, Genentech, Gates Foundation

More information

http://hcbravo.org
@hcorrada